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Abstract 

Self-consistent solutions to a generalized Su-ScliriefFer-Heeger model on a 2-dimensional 
square lattice are investigated. Away from half- filling, spatially inhomogeneous 
phases are found. Those phases may have topological structures on the flux or- 
der, large unit cell bond order, localized bipolarons, or they are simply short-range 
ordered and glassy. They have an universal feature of always possessing a gap at 
the Fermi level. 
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The Su-Schrieffer-Heeger (SSH) model [1] was motivated by the quasi- 1-dimensional 
conducting polymer polyacetylene. The model exhibits Peierls instability at 
the characteristic wavevector 2k where /c^ is the Fermi wavevector that de- 
pends on the band filling. As a result the translational symmetry of the lattice 
is spontaneously broken and a gap is opened up at the Fermi level of the elec- 
tronic spectrum. Since the lattice distortion in the ground state is multiply 
degenerated, topological mid-gap soliton states can be formed at the domain 
boundaries. In-gap polaron states can also be formed as soliton-antisoliton 
bound states. At the empty band limit, bipolaron states are favored. The 
above mentioned states as the elementary excitations of the system and their 
experimental consequences were reviewed in detail in an article by Heeger 
et al. [1] In dimensions higher than one, for instance in a 2-leg ladder, a new 
mid-gap soliton named the "twiston" was discovered [5]. The ground state in 
the half-filled 2-dimensional (2D) square lattice was also studied within an as- 
sumption of a small unit cell [2,3,6], or using some more elaborated numerics 
[4]. In general, studies of the model in higher dimensions are relatively fewer. 
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and knowledge of the possible phases is limited. Therefore investigations of 
the phases in some higher dimensional SSH-type models are desirable. In this 
paper, we have extended these studies to a generalized SSH model in two 
dimension (2D). It is found that generally speaking, the types of phase are 
much more diverse in 2D. Some of them appear to be nice generalizations of 
the above-mentioned states in one dimension (ID). In particular, the transla- 
tional symmetry is also spontaneously broken and the Fermi levels are gapped 
in many phases. Thus this generalized SSH model represents an insightful 
generalization of the conventional model and merits serious consideration. 

We start by considering the conventional SSH model in the adiabatic limit. 
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where t is the hopping, a is the electron-lattice coupling, k is the interatomic 
elastic constant, and 0ji = a — |ri — rj| is the deviation of an interatomic 
distance from its equilibrium value a. Every quantity here is real-valued and 
a and k are positive. To find the zero temperature configuration of the lat- 
tice, 0i+n,i's are treated as variational parameters and adjusted (through some 
procedures such as the steepest-descent method) to minimize the total en- 
ergy (Hssu)- Alternatively, one can recast the minimization problem into an 
equivalent self-consistency problem using d{HssH) /d4>ij = 0, which leads to 
0ij = {2a/ k) < Cio-^Cjo- > and 0i+n,i is now a bond parameter to be solved 
by numerical iteration. Defining $ij = a^ij and A = Aa^/n, the problem is 
reformulated as 



-f^SSH = - E E + ^i+n,i)Cia'^Ci+„,^ 

icr n=±x,±y 

+ XE E l^i.i+nr, (2) 

i n=x,y 

^ij = ^ E < Cia^Cj^ > . (3) 

In this form, it is easy to see that the only quantity characterizing the problem 
is t/ A or Kt/{4:a'^), which is dimensionless. 

In the formulation Eq. 2 and 3, it will be meaningful to extend $ij from the 
real number regime into the complex number regime. The physical appeal of 
the extension is that the model now addresses a class of more general problems 
by the following reason. Generally, the parameter $ij represents an overlap in- 
tegral < > describing off-site electronic correlations, with H' and 
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Fig. 1. The 2-sublattice ansatz. The 2D square lattice is patched up by centering 
this cross at sites belong to one of the two sublattices. The complex-valued nature 
is indicated by an arrow, which means that ^^^^^ = $^f^ is the reversed arrow and 
may not be identical to ^^f^- 

IV'i/j > being the off-site part of the total Hamiltonian and the local atomic 
wavefunctions respectively. Since in addition to the electron-lattice interaction, 
a renormalized low-energy effective theory may contain other many-particle 
effects manifested in the overlap integral, there is no reason for $ij to be nec- 
essarily real when these effects are included. This extension thus goes beyond 
the electron-lattice interaction and it may be viewed as a prototype of this 
class of problems. Indeed, a four-fermion interaction is generated at integrat- 
ing out $ij, and therefore the self-consistency equation that $ij respects can 
also be regarded as a mean field equation. In fact, the model is formally the 
same as a mean field theory in the high-Tc cuprates [7] . In view of this connec- 
tion, we have in this paper investigated this generalized SSH model in detail. 
On one hand, our study is a generalization of the conventional SSH model in 
2D; on the other hand, it is also a study complementary to those mean field 
studies of the high-Tc system, where a different parameter regime (larger t/\) 
was emphasized. 

It is noted that the generalized hamiltonian bears the same form as a mean 
field hamiltonian in the high-T^ cuprates [7], and presumably they should 
bear the same solution. In the high-Tc context, a 2-sublattice (2SL) ansatz 
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(see Fig. 1) designed to capture any (tt, tt) -instability was adopted in Ref. [7] 
for the solution. Such a high translational symmetry solution might be appro- 
priate in that context, since other effects such as fluctuation might at work 
and preserve the symmetry to a high degree, but in our case, where we expect 
soliton or polaron solutions, such an ansatz at the outset is clearly inappropri- 
ate. Nevertheless, we still have reproduced the solution within that framework 
and compare with the result of our approach which allows spatially nonho- 
mogeneous solutions. The 2SL ansatz assumes the distribution of $'s on the 
lattice has a repeating unit as shown in Fig. 1. Different types of solution 
were found [7], they were the staggered- Peierls (SP) with complex ^^^^'s and 
$2SL ^ ^2SL ^ ^2SL > Uniform staggered-flux (u-SF) with 

all the ^^^'"'s are equal and have imaginary components; uniform real-bond 
(u-R) with all the $2^^'s are equal and real; kite (K) with all the ^^^'"'s are 
real and either $^sl ^ ^2SL _^ ^2SL ^ ^2SL ^2SL ^ ^2SL _^ ^2SL ^ ^2Sl_ 

We have found the zero temperature self-consistent solutions iteratively on 
periodic boundary finite square lattices. Each bond is treated as an inde- 
pendent parameter. Meanwhile, phases and energies within the 2SL ansatz 
are also solved and compared with the result of our unrestricted search. Our 
unrestricted-bond iteration converges to the 2SL form only at two regimes, 
(i) X ~ and t > ~ O.IA, where the uniform-SF is located at; and (ii) 
large t and x, where the uniform real-bond phase is located at. x is the hole 
concentration between and 1. 

Away from the above regimes, our solutions generally do not fit in a 2SL 
ansatz. Some representative phases which we will discuss and their energies 
(per site) Eg are grouped here beforehand. In the brackets, energies E^f^ and 
phases from the 2SL ansatz are also given for comparison. We obtain dimer-box 
(D-B) glass phase at, t = 0, a; = 0.3, with Eg = -0.1750A [E^^^ = -0.1608A, 
K]; t = 0, X = 0.5, with Eg = -0.1250A [E^^^ = -0.1077A, u-R]. We find 
striped-staggered-flux (s-SF) at, t/X = 0.2, x = 0.1, with Eg = — 0.5017A 
[E^SL ^ _o.4991A, u-SF]; t/\ = 0.2, x = 0.2, with Eg = -0.4784A [E^^^ = 
-0.4733A, u-R]; t/\ = 0.4, a; = 0.1, with Eg = -0.8134A [E^^^ = -0.8088A, 
u-SF]. Wigner lattice of bipolarons (WBP) is found at, t/X = 0.1, x = 0.75, 
with Eg - -0.1267A [E^^^ = -0.1232A, u-R]. A bond order wave (BOW) 
which wc have named Real-staggered-box (RSB) is found at, t/X — 0.1, x — 
0.5, with Eg = -0.2402A [E^^^ = -0.2390A, u-R]; t/X = 0.05, x = 0.5 with 
Eg = -0.1798A [E^i^^ = -0.1734A, u-R]. Solutions of lower translational sym- 
metry, when they exist, are found to have lower energies and the percentage 
of difference is more prominent at small t but intermediate x. Typical total 
density of states (DOS) of these phases are shown in Fig. 2. 

Dimer-box glass: On the t = axis, a solution of high degeneracy can be 
found analytically. It is straightforward to verify that a lattice arbitrarily 
filled with disjointed dimers and boxes is a solution [13]. A dimer [7,12] is 
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Fig. 2. Total density of states at some representative dopings and hoppings (A = 1): 
The shaded are occupied levels and the horizontal axis of each box is the energy 
from -2.5 to +2.5. 

a doubly filled nearest-neighbor bond [19] with |$| = A/2. A box [14] is a 
plaquette with bonds 
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and a real 



$12*23*34*41 < 0, 



(6) 



where 1, 2, 3, and 4 are consecutive corners of the plaquette. The lowest two 
levels of a box are filled with 4 electrons. Both dimer and box have an average 
of 1 electron per site, energy —A/4 per site, and only two levels at ±A/2. At 
X — 0, the lattice is fully filled with dimers and boxes. At x > 0, a solution 
can be simply obtained by plucking off dimers or boxes [15] from the lattice 
and the energy will be just — (1 — a;)A/4. Creating empty sites also creates zero 
energy states, therefore the spectrum of a doped zero-t system is a three-level 
structured DOS D{£) = x6{£) + (1 - x)[6{e + A/2) + 6{£ - A/2)]. Since the 
arrangement of the dimers and boxes is arbitrary, this phase is in general a 
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glass of bond and plaquette centered electron charges. But it is important to 
note that due to the degeneracy condition Eq. 5, a box can be continuously 
deformed into two dimers or vice versa without any energy barrier intervening. 
Therefore this glass might not be a stable phase when fluctuation effect is taken 
into account, and some translational symmetry is expected to be restored. 

It can be shown analytically that the D-B glass has a lower energy than the 
u-R phase at t = and a; — 1. At a; — > 1, the kinetic energy of a u-R 
phase is simply obtained by filling in the band bottom, i.e., —2D{t + $)(! — 
x), and the potential energy is 2£)$^/A, where D — 2 is the dimensionality. 
Minimizing the total energy gives $ = {1 — x)X/2, hence energy 
-2tD{l - x) - {XD/2){1 - xf. It is of - x)^] at t = 0, while the D-B 
glass has Eq = —(1 — a;)A/4 of — x)]. Comparing E^f^ and Eq also 
gives a scale t — X/ {8D) at which the D-B glass is destabihzed and there is a 
crossover from short-range-order (SRO) to long-range-order (LRO). 

At intermediate x, one can only obtain E'lf^ numerically. We have studied, 
e.g., X = 0.1, 0.3, and 0.5, and have found that they have energies (which are 
given before, and are stable at lattices 40x40, 80x80, and 120x120) substan- 
tially higher than the D-B glass. A feature to note is that the 2SL phases are 
ungapped, whereas the three-level D-B glass is obviously gapped. The levels 
expand into bands at small but nonzero t and the gapping remains, as a vestige 
of the zero-t phase. 

Striped-staggered- flux: At the regime < x < 0.2 and ic(~ O.IA) < t < 
0.6A, the holes are found to arrange into stripes along the lattice axes [16,17] 
(say y-axis, see Fig. 3). Whereas uniform-SF is obtained within the 2SL ansatz. 
The stripes at the meantime are also topological domain walls separating an- 
tiphase SF domains [17]. Two types of domain walls are found, one has a 
column of neutral-flux plaquettes inserted, the other has a column of plaque- 
ttes deleted. It is remarkable that the electron depletion can be as large as 
~0.2 and has a short coherence length of only a few lattice constants. We also 
make an observation that it has one hole per unit length [18]. The spectrum 
[see Fig. 2(b), (e), and (f)] is nodally gapped at the Fermi level and the gap 
maxima at (0, ±7r) and (±7r, 0) [7,12,8,9,10,11] manifest themselves as the van 
Hove singularity peaks in the DOS. Local DOS shows that the mid-band has 
a great part composed of states localized along the stripes (see Fig. 3). 

It is interesting to check if such striped-SF phase have energies lower than 
that of the 2SL phases. For some commensurate fillings like x — 0.1,0.2, 
and 0.3, we have obtained straight hole-stripes along one of the lattice axes. 

In such cases, the fact that the striped-SF phase have lower energy can be 
justified in the thermodynamic limit as follows. Owing to the translational 
invariance along one axis, we may equivalently seek for such solution within 
an ansatz with translational symmetry = ^i+2y,i+a+2y,C( — x,y [20]. 
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Fig. 3. Topological stripe: (Top) The bond current (4f/A)Im$ij, orbital moment 
(average of directed bond currents round a plaquette) , and (middle) electron density 
profiles across an ID antiphase domain wall {t/\ = 0.2, x = 0.1). Bond thicknesses 
are proportional to the current sizes and circle areas to the moment sizes. Maximum 
bond current here is 0.106/s, and maximum moment is 0.105/s. The domain wall 
shown here contains neutral plaquettes. (Bottom) Typical local DOS at hole-rich 
and hole-poor sites are also shown. 

Due to the reduction of numerical effort (diagonalizing a hamiltonian on a 
NxXNy lattice is now separately diagonalizing Ny/2 hamiltonians on NxX2 
lattices), self-consistent solutions can now be pursued on much larger lattices. 
We have justified that the striped-SF at, e.g., t/X = 0.2, x = 0.1,0.2, and 
t/X — 0.4, X = 0.1 have energies (which are given before) lower than that of 
the 2SL solutions. The energies are stable at lattices from 40x40 to 120x120. 



Note that the stripes here are not a result of frustrated phase separation 
since we assume no repulsion between holes. They are more appropriately 
understood as soliton states created to accommodate the holes, analogous to 
the soliton states in the ID SSH model [1] . By segregating out the holes, the 
SF order has survived into higher hole concentrations where it would not has 
existed within the 2SL ansatz. 



LRO-SRO crossover: The nonzero but small t regime, t < tc ^ O.IA is a 
crossover region between the large-t regime of LRO and the zero-i regime 
of degenerated glass. SRO dominates over this regime and there are lots of 
almost-degenerate local minimum solution. Phases in this regime are generally 
glassy, have complex bonds at x < 0.3, and real bonds at x > 0.3. Charge 
distribution is correlated with the local bond order and also nonuniform in 
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Fig. 4. (Left) Real-staggered-box BOW with uniform charge distribution at 

t/\ = {).l,x = 0.5. (Right) Wigner lattice of bipolarons at t/X = 0.1,, x = 0.75. 
Bond thicknesses and circle areas are drawn proportional to the bond magnitudes 
and electron densities respectively. 

general. Basically, regions that are crowded with electrons may have complex 
bonds and regions that are not, have essentially real bonds. In our calculation, 
such glassy states found by iteration do have energies lower than that of the 
2SL solution and their spectra are always gapped at the Fermi level. 

At a; ~ 1 of this crossover regime, the obtained state contains uncorrelatedly 
scattered bipolarons. A bipolaron is a localized state occupied by two electrons, 
with nonzero real bonds at its vicinity. Decreasing t shrinks them into dimers, 
or increasing electron density turns them into an almost-hexagonal closed- 
packed Wigner lattice (WBP) (see Fig. 4). In the well-isolated bipolarons 
regime, the energy can be readily obtained as the sum of the energies of the 
individual bipolarons. For the Wigner lattice, for instance at t/A = 0.1, a; = 
0.75, we obtain the same energy (which is given before) at lattices from 20x20 
to 32x32, and it is lower than the energy of the 2SL solution (also given 
before) . 

Two exceptions out of these glassy phases are the BOWs with uniform charge 
distribution named "real-staggered-box" (RSB) (see Fig. 4) found aX x — 0.5, 
and "box-staggered-flux" (B-SF) at a; = 0. RSB found at t/X = 0.1 and 0.05 
have energies (which are given before, and are stable at lattices from 40x40 
to 120x120) lower than that of the u-R. B-SF is a phase with SF order and 
stronger bonds distributed like a square array of boxes (resembles those boxes 
in Ref. [7]), and it also has energies lower than that of SP. 
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Fig. 5. Approximate phase diagram: B-SF: Box-staggered- flux (gapped). u-SF: 
Uniform staggered-flux (ungapped when x > 0). s-SF: Striped-staggered-flux 
(gapped). u-R: Uniform real-bond (ungapped). D-B: Dimer-Box glass (gapped). 
The dimer and box order are exact only at t = 0. g-R: Phases with glassy real-bond 
and charge density (gapped), and are not characterizable by a simple order. RSB: 
Rcal-staggered-box BOW with uniform charge density (gapped). WBP: Wigner lat- 
tice of bipolarons (gapped). i-BP: Isolated and uncorrelated bipolarons (gapped). 

When X < 0.3, bonds are complex but remain glassy. The spectrum has a 
narrow in-gap band consists of almost-localized states. These in-gap states are 
fundamentally different from those of the striped-SF at larger t. The latter are 
soliton states supported by the LRO, while the former are states localized in 
the bond glass. 

In summary, an approximate phase diagram is given in Fig. 5. Spatially non- 
homogeneous phases are shown to be energetically more favorable in the gen- 
eralized SSH model in 2D away from half-filling, and the electronic spectrum is 
believed to show a similar gapped Fermi level as that in the ID conventional 
SSH model. The doped holes also go into the in-gap states and topological 
object may also be formed, but the 2D phases are much more complicated 
than the simple polymerizations in ID. Those phases may be superstructurcd 
or short-range-ordered. At large t/\ and large x, the translational symmetry 
is preserved and this is fundamentally different from the ID case, where the 
Peierls instability always in effect. 
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